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Abstract The cause of enhanced acoustic power surrounding active regions, the 
acoustic halo, is not as yet understood. We explore the properties of the enhanced 
acoustic power observed near disk center from 21 to 27 January 2002, including 
AR 9787. We find that (i) there exists a strong correlation of the enhanced high 
frequency power with magnetic-field inclination, with greater power in more 
horizontal fields, (ii) the frequency of the maximum enhancement increases along 
with magnetic field strength, and (iii) the oscillations contributing to the halos 
show modal ridges which are shifted to higher wavenumber at constant frequency 
in comparison to the ridges of modes in the quiet-Sun. 

Keywords: Helioseismology, Observations; Active Regions, Magnetic Fields; 
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I. Introduction 

The application of local helioseimology to probe the nature of solar active regions 
is critically important in understanding their origin and evolution. The influence 
of magnetic fields on incident acoustic waves is complex, but is coming under in- 
creased scrutiny due to advances in numerical simulations of wave propagation in 
magnetised plasmas (e.g.lKhomenk o and Co llados, 2006 ; Shelyag, Erdelyi ,"and Thompson[ | 



l2006[[Parchevsky and Kosovichev|l2007j|Cameron, Gizon, and"T3 uvall, 2008; Rempel, Schiissle r, and Knolker 



2009}. A successful model should reproduce all of the observed features, including 
changes in the surface amplitudes of waves in the vicinity of active regions. 

The acoustic halo is an observed enhancement of high frequency acoustic 
power surrounding regions of strong magnetic field. This enhancement was dis- 
covered independently by IBraun ~eTaL\ (|1992aj) , [B rown et aTl ([1992) and 
IToner and Labontel ([T993 ) . Observations made using Ca II K-line images showed 
these enhancements extending tens of Mm beyond active regions (Br aun et a/., 1992a| | 



Toner and Labonte, 1993} . IBrown et al . (1992) showed patchier and more con- 



fined halos using Doppler velocity measurements taken in the Fe I A5576 line, and 
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with higher spatial resolution. They also showed a peak in the power at 5.5 mHz. 
IHindman and Brownl (|1998|) and later [Jain and Haber (2002) demonstrated 
that halos seen in SOHO/MDI (Scherrer et a/., 1995ft Doppler grams tend to be 
prominent in intermediate magnetic field strengths of 50-250 G and are absent 
in continuum intensity observations. Acoustic halos have been shown to be more 
prominent in Ca II K-line power maps and less visible in the Doppler Ni I line 
shifts ( Ladenkov et a/., 2002 ). IFinsterle et al.\ (|2004p found that the acoustic 
halo increases in size with increasing height into the canopy over active regions. 



On the other hand, Muglach (2003]) found that high-frequency oscillations in 



TRACE UV filtergrams show deficits in magnetic fields, rather than enchance- 
ments. There is no established explanation for the enhanced power, although 
several ideas have been suggested (see Section Hj). 

In this paper we present an analysis revealing several new properties of acous- 
tic halos and of surface acoustic power in general. We use data obtained by 
SOHO/MDI over 7 days containing several large active regions. In Section [2] 
we demonstrate (i) a strong association of the enhanced high frequency power 
with nearly horizontal field, and (ii) an increase in the temporal frequency of 
the maximum enhancement with field strength. In Section [3] we show that (hi) 
the oscillations contributing to the halos show modal ridges which are shifted to 
higher wavenumber at constant frequency than the ridges due to waves in the 
quiet-Sun. In the Discussion we list possible theories put forward, to date, that 
attempt to explain the existence of the acoustic halos, noting that none of them 
have been shown to explain all of the observed properties successfully. 



2. Correlations between the acoustic power enhancement and the 
magnetic field 

We use 7 days of observations from 21-27 January 2002 covering a significant 
number of active regions including AR 9787 which has already been subjected 
to a comprehensive helioseismic analysis ( |Gizon et a/., 2009 ). We track the mag- 



netograms and Dopplergrams, remapped into Postel's projections, in 24 hour 
blocks centered at the disk center for the middle of that day. The scale of the 
map is O.OO2i?0 « 1.4 Mm per pixel over 700 pixels, so the extent of the maps 
are about 975 Mm in each direction. The acoustic power maps are normalised 
to the quiet-Sun power by removing a spatial 2D 4th degree polynomial fit to 
the quiet-Sun power for each day and frequency band. This fit removes several 
effects including center-to-limb variations in the observed mode amplitudes and 
spatial fluctuations in the instrumental MTF. The latter includes variations of 
focus across the CCD. 

The quiet-Sun is defined to be those pixels where the total magnetic field 



strength, \B\ = y B% + B^ + B%, is less than 30 G. Standard (potential field) 

methods were applied to reconstruct the vector magnetic field, (B x , B y , B z ), at 
the Sun's surface from the line-of-sight components given by the SOHO/MDI 
observations (see ISakurail (|1982|) for the more general problem of extrapolat- 
ing the field into the overlying corona). The ^-direction is defined as positive 
in the direction of solar East perpendicular to the meridian, the ^/-direction 
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is defined as positive towards solar South perpendicular to the equator and 
the radial magnetic field is represented by B z . Keeping in mind that poten- 
tial field extrapolations are not the best representation of the actual magnetic 
field (Wiegelm ann et a/., 200 5), we avoid using the strong fields (for example, in 
sunspots) where such extrapolations are known to be particularly unreliable. 

The first four panels of Figure [1] show acoustic power in 1 mHz bands centered 
on (a) 3 mHz, (b) 4 mHz, (c) 5 mHz and (d) 6 mHz for the 24 January 2002. 
The sunspot in AR 9787 is seen close to the center. The bottom two panels show 
intensity continuum (e) and the line-of-sight magnetic field (f). In the lower 
frequency bands, (a) and (b), we see suppression of acoustic power coincident 
with the location of magnetic field. In (c) we begin to see an enhancement of 
power surrounding the suppression in plage regions, and in (d) we see predom- 
inantly enhanced power except in the very strongest of magnetic fields, within 
the sunspot itself. It is clear that both the suppression (at lower frequencies) and 
enhancement (at higher frequencies) are spatially correlated with the magnetic 
field, even in the smallest plage regions. It is also interesting to note from panels 
(d) and (f) that the strongest 6 mHz enhancements are located between regions of 
opposite polarity (indicated by the arrows), leading us to explore the dependence 
upon magnetic field inclination. 

Figure [2] shows the vector magnetic field, (a) B X1 (b) B y and (c) \B\ derived 
from the line-of-sight observations in Figure [T] (f). From these maps we made 
scatter plots of the acoustic power against the field inclination for each of the 
frequency bands shown in Figure[TJ The field inclination, 7, is defined as tan(7) = 



B z /Bh where Bh = y B% + By is the horizontal field strength and B z is the field 
parallel to the radial vector. Thus, purely horizontal field is defined when 7 = 0°, 
7 < 0° when B z < and 7 > 0° when B z > 0. Figure [3] shows the results for 
(intermediate) magnetic field strength, 150 < \B\ < 350 G, over each of the four 
frequency bands centered at 3, 4, 5, and 6 mHz. For all frequency bands, there 
is a pronounced peak in power at 7 = 0°. This is true for the lower frequencies, 
at which the power in the magnetic regions shows suppressed levels compared to 
the quiet-Sun, but is more pronounced at 5 mHz and 6 mHz. Acoustic power in 
these intermediate-strength magnetic pixels at 5 mHz can be either suppressed 
(in more vertical fields) or enhanced (in more horizontal fields), compared to the 
quiet-Sun. At 6 mHz all of these pixels show enhanced power. 

Figure H] shows the variation of acoustic power in the 6 mHz bandpass with 
inclination over a wider range of total magnetic field strength. For clarity, only 
the acoustic power, averaged over bins in 7 comprising 1% of the total number 
of points, is shown for each magnetic field strength range. It is clear that the 
trend of increased power in horizontal fields is present for all ranges of field- 
strength shown, although the results become noisier at higher strengths due to 
poorer statistics. The greatest enhancement (50% above the quiet-Sun power) 
occurs for field strength between 300 and 400 G, and falls off with either greater 
or lesser strengths. There is a power suppression at the highest field strengths, 
occurring mostly in the sunspots. 

To visualise the dependence of the normalised power on both the field strength 
and inclination, we compute the average power in 10 G bins of field strength 
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Figure 1. 24 January 2002 AR 9787 acoustic power maps at (a) 3 mHz, (b) 4 mHz, (c) 
5 mHz, (d) 6 mHz, with the intensity continuum (e) and the line-of-sight magnetic field (f). 
The acoustic power maps are normalised to one in the quiet-Sun. The arrows in panel (d) 
indicate bright halos located between opposite polarity regions as seen in panel (f). The length 
of each side of the shown observations are « 975 Mm. 



and 10° bins of inclination. Figure [5] shows the 2D distribution of the mean 
normalised power for 1 mHz bands centered on (a) 3 mHz, (b) 4 mHz, (c) 
5 mHz and (d) 6 mHz. The distribution of enhanced power is symmetric about 
the horizontal field at 5 and 6 mHz, with most of the enhanced power occurring 
within a narrow range of inclination and field strength. 
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Figure 2. 24 January 2002 AR9787 potential field extrapolation (a) B x , (b) B y , (c) \B\, from 
to 1000 G and (d) 7. B z is not shown here, but it is the same as the line-of-sight component 
in Figure ^(f). The length of each side of the shown observations is ~ 975 Mm. 

We explore the possibility that the functional dependence shown in Figure [5] is 
sufficient to predict the appearance of the acoustic halo using only the extrapo- 
lated magnetic field components. Figure [6] (left column) shows the reconstructed 
acoustic power maps using the magnetic field for 24 January 2002 (as shown in 
Figure [2]), and the 2D distributions illustrated in Figure [5] The right column of 
Figure [6] shows the actual normalised acoustic power map for 24 January 2002. 
Each row of Figure [6] represents a different frequency band. The reconstruction 
is qualitatively successful with the correlation coefficient, C, between the maps 
given on the right-hand-side. The correlation is best for low frequencies. Recall 
that these maps are created only for the magnetic field map of one day (24 
January 2002) but using Figure [5] which is an average of power over all seven of 
the observed days. Even so, the correlation coefficient for using power from the 
same day is very similar (C is 0.74, 0.73, 0.55, 0.53 for 3, 4, 5 and 6 mHz bands 
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Figure 3. Normalised acoustic power against inclination, 7, for magnetic field strength 
150 < \B\ < 350 G for (a) 3 mHz, (b) 4 mHz, (c) 5 mHz and (d) 6 mHz. The blue curve 
is the binned average of 200 points. The vertical line indicates horizontal field. 



respectively). It therefore appears that the acoustic halo is largely a property of 
the local surface magnetic field. 

There is some indication from Figure [5] that the peak in the mean normalised 
power at 6 mHz occurs at a higher field strength than the peak in the distribution 
at 5 mHz. To explore this variation further, we construct a new set of normalised 
acoustic power maps with narrower frequency bandpasses. From the symmetric 
properties of the acoustic power around 7 = 0° we can comfortably work with 
|7|. Figure [3 shows the mean acoustic power in 10 G bins of magnetic field and 
0.1 mHz bins of frequency for 0° < I7I < 30° (left) and 30° < I7I < 90° (right). 
A clear trend of increasing peak frequency with \B\ is evident in both ranges of 
inclination. The well known suppression of power at 3-4 mHz, which is stronger 
at greater field strength, also shows up prominently in Figure [71 
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Figure 4. Mean 5.5 - 6.5 mHz acoustic power against inclination in ranges of magnetic field 
strength as indicated. The points have been averaged in bins of inclination comprising 1% of the 
total number of points in each field-strength range and the standard deviation is over-plotted. 
Each curve, except the top, has been shifted down for ease of display, where the dashed curves 
represent a value of 1 for each magnetic field band. 



3. Wavenumber dependence 

We now turn to examining the nature of the enhanced power oscillations in 
terms of changes in the frequency and wavenumber domain. To isolate the 
power spectra of the halo regions we apply spatial masks to the observations 
before performing a 3D Fourier transform. Since the power spectra are quite 
sensitive to the window function, or mask, it is necessary to compare the power 
spectra of observations of the halos with those of quiet-Sun using identical masks. 
Because of the large filling factor of the active regions in the 2002 January 
data this comparison nominally requires an additional 7 days of observations of 
predominantly quiet-Sun. This entails using additional MDI observations from 
a substantially quiet epoch. For the quiet-Sun comparisons used here, we chose 
observations (hereafter referred to as the "2003 observations") from 24 to 27 
July 2003 and 18 to 20 October 2003. We refer to the 7 days of data for the 
active period from 21 to 27 January 2002 as the "2002 observations". We track 
and Postel-project the 2003 observations in an identical manner to the 2002 
observations and then pair the 2002 and 2003 observations in days according to 
Table [U 

Unfortunately, focus changes in the MDI instrument were made between the 
2002 and 2003 observations. In addition, the time span between the two datasets 
is substantial enough that possible solar cycle changes in mode properties could 
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Figure 5. Mean acoustic power in 10° bins of inclination, 7, and in 10 G bins of magnetic 
field strength in 1 mHz band passes centred at (a) 3 mHz, (b) 4 mHz, (c) 5 mHz and (d) 
6 mHz averaged over 7 days (21 - 27 January 2002) where 7 = 0° is horizontal field. 

confuse our power spectra comparisons. Thus it is necessary to look for, and rule 
out, possible changes in the power spectra of the quiet- Sun over the same two 
epochs which might be due to either instrumental or solar causes. We therefore 
create two sets of masks; the first of which (called the "halo masks") is intended 
to isolate the halo regions in the 2002 observations, and the second ("quiet- 
Sun masks") is intended to isolate the quiet regions in the 2002 observations. 
An additional complication is that the 2003 observations also include magnetic 
regions which must be eliminated from the comparison. To do this, the halo 
masks for each pair of days include only pixels with 100 < \B\ < 350 G and 
|7| < 20° in the 2002 observations and with \B\ < 30 G in the 2003 observations. 
The set of quiet- Sun masks is constructed for each pair of days from pixels with 
I B I < 30 G in both the 2002 and 2003 observations. This excludes any magnetic 
pixels present in the 2003 observations, and ensures that the masks used for both 
sets of days are identical. 
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Figure 6. The reconstructed acoustic power maps for 24 January 2002 are shown in the left 
column and the right column shows the observed acoustic power map. Each frequency band 
(3 mHz, 4 mHz, 5 mHz and 6 mHz) is shown from top to bottom row. The reconstruction 
is qualitatively successful with the correlation coefficient, C, between the original and recon- 
structed power maps given on the right-hand-side. The length of each side of each map is 
& 975 Mm. 
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Figure 7. Mean acoustic power in 0.1 mHz bins of frequency and in 10 G bins of magnetic 
field strength averaged over 7 days (21 - 27 January 2002) for 0° < \j\ < 30° (left) and 
30° < |7| < 90° (right). The lowest power contour is lighter in colour for visibility. 



Table 1. The dates of observations for pairing the 2002 observa- 
tions (January) and the 2003 observations (July and October). 
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Using these masks we compute the azimuthally summed power spectrum for 
each pair of days. The 7-day average of the 2002 and 2003 spectra are shown in 
Figure [U It is apparent that the power in the 2002 spectra of the halo regions 
is greater, and the ridges are more defined, than the spectra in 2003. Moreover, 
there appears to be a shift towards higher k in the mode ridges in the halo power 
spectra which is particularly visible at higher temporal frequency. 

To show this more clearly, we plot cuts through the spectra in Figure [9] at 
6.3 mHz (panel a) and at 5.9 mHz (panel b), as indicated by the vertical lines 
in the power spectra (Figure [8]). To facilitate the comparison at 5.9 mHz the 
quiet- Sun cut has been multiplied by a factor of 1.21. From these plots it is clear 
that there is a shift of the ridges in the 2002 halo-masked power spectrum to 
higher k. 

Similar shifts of the ridges between the 2002 and 2003 quiet-Sun regions, if 
present, would indicate some artifact or other effect which is not unique to the 
halos. Figure [10] shows cuts through the power spectra obtained with the quiet- 
Sun masks. It is apparent that some differences in the magnitude and distribution 
of power, but not frequency shifts, between the two epochs are observed. The 
former is not unexpected given the known focus changes between the two epochs. 
However, it is clear that we do not see a shift in the position of the ridge peaks 
at constant frequency as is seen in Figure [9] This strengthens the likelihood of 
real differences in the modal properties between halo and quiet regions. 
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Figure 8. The azimuthally summed power spectra (in arbitrary units) for the halo masks 
applied to the 2002 (active region) observations and averaged over all days (right panel) and 
for the same masks applied to the 2003 (quiet-Sun) observations averaged over all days (left 
panel). The vertical lines correspond to the cuts plotted in Figure [9] 

4. Discussion 

We find properties of the acoustic halo that confirm prior analyses, but also dis- 
cover several heretofore unknown properties. The halos are locations of acoustic 
power increase up to 140% of the quiet- Sun value at higher frequencies and 
are strongest at intermediate magnetic field strength (150 to 350 G). These 
results are in general agreement with previous studies QJain and Haber, 2002 ; 
Hindman and Brown, 1998| . We have discovered a clear association with near 
horizontal (±30°) magnetic field. The dependence of acoustic amplitude on 
field inclination occurs over frequencies between 2.5 and 6.5 mHz, and for field 
strengths between 100 and 800 gauss, but is particularly strong between 5 and 
6.5 mHz and at the aforementioned intermediate field strengths. In addition, we 
show that the frequency of the peak enhancement depends on the field strength, 
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Figure 9. Cuts through the power spectra (in arbitrary units) shown in Figure [8] at constant 
frequency. The top panel (a) shows a cut through the power spectrum at 6.3 mHz for the 
halo-masks applied to the 2002 (active region) observations (blue) and a cut through the the 
2003 (quiet-Sun) power spectrum (black). The same is shown in the bottom panel (b) except 
at a frequency of 5.9 mHz. At 5.9 mHz, the 2003 (quiet-Sun) power spectrum has been boosted 
by 1.21 to allow an easier comparison with the 2002 spectrum. Some smoothing of the curves 
was performed for the purpose of the plot. 



with the peak frequency increasing slightly with field strength, and that in these 
regions of enhanced power the ridges of the acoustic power spectrum are shifted 
towards higher k. Qualitatively there appears to be a larger shift for higher k 
indicating that this is a shallow surface effect. For constant k the phase speed of 
the ridge is reduced, indicating an increase in travel times, that could be caused 
by either a slower wave speed or a longer path distance. 
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Figure 10. The same as Figure [9] except the cuts are taken through the power spectra after 
applying the quiet-Sun mask set. The top panel (a) shows smoothed power cuts at 6.3 mHz 
for quiet-Sun regions in the 2002 observations (blue) and quiet-Sun regions from the 2003 
observations (black). The same is shown in the bottom panel (b) except at a frequency of 
5.9 mHz. 



Any theory attempting to explain the acoustic halo or fully explain the inter- 
acation of waves with the magnetic field needs to incorporate these properties. 
A quantitative study of changes in the power spectrum requires a much more in 
depth analysis than our rudimentary study discussed here since, amongst other 
things, we have made no attempt to remove or account for distortions due to 
foreshortening and projection effects from the power spectra. On the other hand, 
we are confident through the multiple comparisons performed here with different 
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masks and solar regions that the ridge shifts are qualitatively real and are clearly 
associated with the halo regions. 

Various authors have pointed to a few different mechanisms to explain the 
halo. Some of the original discovery papers suggested enhanced acoustic emis- 
sion QBraun et a/., 1992b} |Brown et a/., 19 92). More recently, IJacoutot et al.\ 
(|2QQ8p have studied the effects of vertical magnetic fields on solar convec- 
tion. The convection occurs on smaller scales and higher frequencies than in 
quiet- Sun regions, and shifts the oscillation power to higher frequencies. It is 
not as yet understood whether these properties are dependent upon the in- 
clination of the field. However, the lack of halos in continuum observations 
(Hindman and Brown, 1998) and the relative lack of association with emission 
sites, as identified with seismic holography analyses sensitive to propagating 
waves flDonea, Lindsey, and Braun, 2QQQ[ ) casts doubt on this interpretation. 



Hanasoge (2008) claims that halos may be more likely related to scattering 



effects. Roughly in agreement with the findings in this paper, he also sees the 
halo as being confined to magnetic field inclinations \j\ < 40° and an increase 
of energy in the higher wave number regime. IKhomenko and Colladosl (|2009|) 
propose that high-frequency fast mode waves refracted in the higher atmosphere 
due to a rapid increase in the Alfven speed are responsible for the seismic halos. 
They get a 40-50% increase in power of high frequency vertical velocity in regions 
of intermediate field strength, although, like [Hanasoge (2008]) , they find that 
the horizontal velocity component of the halo is stronger. However, it is difficult 
to completely reconcile the observations with these simulations, since the latter 
show only a weak halo being generated in the vertical velocity component (up to 
a maximum of 5%) and the strongest enhancement (up to 35%) being generated 
in the horizontal velocity component. If the halo were strongest in horizontal 
velocities there would be a line-of-sight dependence with the halo power being 
stronger at larger heliocentric distance. This has not as yet been observed. 

IKuridze et al . (2008) suggests this enhanced power is caused by high fre- 
quency waves trapped beneath the bipolar magnetic canopy. The low frequency 
waves are reflected below the photosphere and so are not trapped. In regions of 
open field, the high frequency waves are permitted to propagate upwards and 
are lost to the atmosphere. These field free cavities have granular dimensions 
(« 0.5 Mm), but with SOHO/MDI we must be observing many of these cavities 
within the resolution which is three times this size, and these cavities must be 
widely distributed to cover the area of the halo regions. A somewhat similar idea 
is put forth by |Muglach, Hofmann, and STa ude (2005]) . They suggest that the 
lack of haloes in the TRACE observations ( Muglach, 2003[ ) caused by reflections 
(and mode conversion) due to an overlying inclined magnetic field. 

ISimoniello et al.\ (|2009p have studied Sun-as-a-star oscillations and observe 
enhanced power in high frequency (5.7-6.3 mHz) bands, up to 18%, at solar 
maximum. We believe that this is due to a large surface of the Sun being 
covered by plage and the associated enhanced power, as presented in this paper, 
associated with the increased magnetic activity. 

Understanding the properties of surface acoustic power is important in mod- 
eling helioseismic signatures due to subsurface perturbations. It has also been 
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suggested that acoustic power be used as a predictive tool to forecast the emer- 
gence of active regions ( |Hartlep et a/., 2QQ9[ ). The results shown here, namely 
that the surface acoustic power is a strong function of the local magnetic field 
properties, may be useful in the removal of surface contributions from acoustic 
power maps and as a part of a correction to local helioseismic methods (e.g. 
|Lindsey and Braun] ()2005|) , |Couvidat and R ajaguru] ( |2007|) ). 



At this point, we cannot confidently suggest a mechanism explaining all the 
properties of the acoustic halo. However, a more detailed analyses of the charac- 
teristics of the acoustic halo will likely help to pin-point the mechanism behind 
the acoustic halos. The recently launched Solar Dynamics Observatory hosts the 
Helioseismic and Magnetic Imager instrument which has a resolution four times 
better than SOHO/MDI and, along side improvements in numerical modeling, 
may be able to help resolve the nature of acoustic halos. 
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